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The idea of exciton condensation in solids was introduced in 1960’s with the analogy to supercon¬ 
ductivity in mind. While exciton supercurrents have been realised only in artificial quantum-well 
structures so far, the application of the concept of excitonic condensation to bulk solids leads to a 
rich spectrum of thermodynamic phases with diverse physical properties. In this review we discuss 
recent developments in the theory of exciton condensation in systems described by Hubbard-type 
models. In particular, we focus on the connections to their various strong-coupling limits that have 
been studied in other contexts, e.g., cold atoms physics. One of our goals is to provide a ’dictionary’ 
which would allow the reader to efficiently combine results obtained in these different fields. 
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FIG. 1. Schematic picture of excitonic condensation in a 
semimetal. Left: The band structure of a semimetal with 
overlapping bands and no inter-band hybridisation. Right: 
The band structure of an excitonic insulator with a gap 
opened by the Weiss field Aec. 


I. INTRODUCTION 

The description of states of matter in terms of spon¬ 
taneously broken symmetries'’^ is one of the most fun¬ 
damental concepts in condensed matter theory. The 
text book examples are crystalline solids (broken trans¬ 
lational symmetry) or magnets (broken spin-rotational 
symmetry). Besides geometrical symmetries quantum- 
mechanical systems possess internal (gauge) symmetries 
associated with conserved charges. The most famous case 
of spontaneously broken internal symmetry is supercon¬ 
ductivity. Another example of broken internal symmetry 
is excitonic condensation. 

Excitonic condensation denotes spontaneous coherence 
between valence and conduction electrons that arises as 
a consequence of electron-electron interaction. As such, 
it is usually represented by orbital or band off-diagonal 
elements of one-particle density matrix, which possess fi- 
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nite expectation values below the transition temperature 
Tc- It was originally proposed^’^ to take place in the 
vicinity of semimetal-semiconductor transition depicted 
in Fig. 1. In an idealised system where the number of par¬ 
ticles in the valence and conduction bands is separately 
conserved, the excitonic condensation breaks this symme¬ 
try and fixes the relative phase of valence and conduction 
electrons. Since inter-band hybridisation or pair- and 
correlated-hopping terms in the electron-electron inter¬ 
action, which mix the valence and conduction bands, in¬ 
evitably exist in real materials the excitonic condensation 
in the above sense is always an approximation. The clos¬ 
est realisation of this idealised exciton condensate (EC) 
can be found in bi-layer systems. In typical bulk struc¬ 
tures, where the conduction-valence charge conservation 
does not hold and thus the relative phase of valence and 
conduction electrons is not arbitrary, the excitonic con¬ 
densation reduces to breaking of lattice symmetries, e.g., 
leading to electronic ferroelectricity®, and/or gives rise 
to exotic magnetic orders characterised by higher-order 
multipoles. The concept of EC allows unified description 
and understanding of these transitions. 

The theory of excitonic condensation was developed 
in 1960’s and 1970’s with weakly correlated semimet¬ 
als and semiconductors in mind. More recently the 
ideas of excitonic condensation were applied to strongly- 
correlated systems described by Hubbard-type models. 
In their strong-coupling limit, these models lead to spin 
or hard-core boson problems that have been studied in 
other contexts. The purpose of this topical review is to 
summarise the recent work on the excitonic condensa¬ 
tion in Hubbard-type models, discuss the corresponding 
phase diagram and describe connections to various other 
models such as the Blume-Emmery-Grifliths®, bosonic 
t-J or bi-layer Heisenberg^ models. In particular, we 
aim at providing a ’dictionary’ for the names of mutu¬ 
ally corresponding phases found in different models. We 
only briefly touch the active field of bi-layer systems in 
Sec. VB and the spectacular transport phenomena ob¬ 
served there and refer the reader to specialised litera¬ 
ture. We completely omit another major direction of the 
exciton-polariton condensation®. 


A. Brief history of the weak-coupling theory of 
excitonic condensation 

In 1961 Mott® considered metal-insulator transition 
in a divalent material (even number of electrons per 
unit cell) associated with continuous opening of a gap, 
e.g., due to application of external pressure. He argued 
that the one-electron picture cannot be correct, that a 
semimetal with small concentration of electrons and holes 
will be unstable if electron-hole interaction is taken into 
account. Knox^ approached the transition from the in¬ 
sulator side and argued that if the excitonic binding en¬ 
ergy overcomes the band gap the system becomes unsta¬ 
ble. These proposals were put into a formal theory by 


Keldysh and Kopaev®, and des Cloizeaux®® who devel¬ 
oped a weak-coupling Hartree-Fock theory of excitonic 
condensation analogous to the BCS theory of supercon¬ 
ductivity. The role of the pairing glue is played by the 
inter-band Coulomb interaction, which favours formation 
of bound electron-hole pairs, excitons. The early theories 
employed so called dominant term approximation, which 
consists in keeping only the density-density intra- and 
inter-band terms of the Coulomb interaction. This ap¬ 
proximation appears well justified for weakly correlated 
metals or semiconductors where the pair glues arise from 
the long-range part of the interaction. It leads to a large 
manifold of degenerate excitonic states. It is due to this 
degeneracy that the small and so far neglected exchange 
and pair-hopping interaction terms may play an impor¬ 
tant role. Halperin and Rice®® showed that these terms 
select one out of the four following states as the lowest 
one: charge-density wave charge-current-density wave, 
spin-density-wave and spin-current-density wave. 

In 1970’s Volkov and collaborators showed in a series 
of articles®^^®'® that excitonic condensation in slightly off- 
stoichiometric systems leads to formation of a uniform 
ferromagnetic state although the normal state does not 
exhibit a magnetic instability. The idea of an excitonic 
ferromagnet was revived at the beginning of 2000’s in the 
context of hexaborides. Sec. V A. 

II. SPINLESS FERMIONS 

Before discussing the two band Hubbard model (7), we 
look at its spinless version 

77efk + H.c. 

{ij) 

+ {yabb\aj -\- Vbadibj^ + H.c. 

(b) 

i i 

which describes electrons of two flavours (a and b) mov¬ 
ing on a lattice and interacting via an on-site interaction. 
Here, aj (o^) is an operator creating (annihilating) a 
fermion of flavour a on lattice site i, nf = aja^ is the cor¬ 
responding local density operator, and analogically for b 
fermions. We will use 4’a{^) and ^/’^(r) to refer to the den¬ 
sity distribution around the lattice sites in orbitals a and 
&, respectively, when necessary. The sum runs over 

the nearest neighbour (nn) bonds. Model (1) is usually 
called extended Falicov-Kimball model (EFKM), where 
’extended’ implies that both ta,tb yf 0. In the original 
Falicov-Kimball model (FKM)®®, 4 = Vab = Vba = 0, 
only the a-electrons can propagate while the ’heavy’ b- 
electrons are immobile.®® 

Three versions of EFKM have been discussed in the 
literature in the context of excitonic condensation: (i) 
FKM with tb = Vab = Vba = 0, (ii) EFKM without 



3 


cross-hopping Vat = Vba = 0 and (iii) EFKM with cross¬ 
hopping Vab, Vba 7 ^ 0. Since the conservation of the 
charge per flavour plays a central role in excitonic con¬ 
densation we compare models (i)-(iii) from that perspec¬ 
tive. Hamiltonian (1) conserves the total charge for any 
choice of parameters. Since the corresponding t/(l) sym¬ 
metry is not broken in any phase considered here, we will 
not mention this symmetry explicitly in the text. The 
FKM (i) conserves the number of the heavy 6 -electron 
on each site, which gives rise to a local C/(l) gauge sym¬ 
metry. Finite tb in (ii) removes the local gauge invari¬ 
ance, but preserves independently the number of a- and 
6 -electrons, which gives rise to a global C/(l) symmetry 
associated with the arbitrariness of the relative phase of 
the a- and 6 -electrons. Finally, finite cross-hopping in 
(iii) removes this U{1) symmetry. An important special 
case of (iii) is a system with a site symmetry that pro¬ 
hibits an on-site a — b hybridisation. As an example, we 
may consider orbitals a and 6 of different parity, which 
implies Vab = —Vba- In this case, the C/(l) symmetry 
is not removed completely but reduced to a discrete Z 2 
symmetry reflecting the invariance of ( 1 ) under the trans¬ 
formation Ui —)■ —tti, Vab Vba- 

On a bipartite lattice there are additional symmetries. 
Models (ii) with tatb > 0 and tatb < 0 can be mapped 
on each other by ai —7 (—l)*ai^^. In the FC phase the 
map turns ferro-EC order (tatb < 0) into an antiferro- 
EC order (tatb > 0). This property will be important for 
understanding the nature of EC phase in FKM (tah = 
0 ), which lies between the ferro- and antiferro- phases. 
There is also a symmetry with respect to the sign of A 
for |ta| = 141 consisting in the exchange 0 ^ 06 ^. 

Model (ii) with symmetric bands ta = tb is identical 
to a single band Hubbard model, where flavours a and 
6 play the role of the spin variable and A stands for an 
external magnetic field. 

Model (ii) with 4 = —tb is identical to the single band 
Hubbard model with attractive interaction. The corre¬ 
sponding mapping is achieved by particle-hole transfor¬ 
mation in one of the bands 6 ^ o 6 j and reversing the sign 
of the interaction U ^ —U. This transformation maps 
an excitonic insulator onto an s-wave superconductor. 
Note that crystal-field A in the excitonic insulator plays 
the role of chemical potential in the superconductor and 
the chemical potential in the excitonic insulator plays the 
role of magnetic field in the superconductor. This shows 
that deviations from half-filling are detrimental for EC in 
the same way magnetic field is for spin-singlet supercon¬ 
ductivity. The mapping between the superconductor and 
excitonic insulator obviously breaks down in the presence 
of external electro-magnetic field. 



FIG. 2. The T = 0 phase diagram of the d = 2 XXZ model 
constructed from data of Ref. 27. The dotted line marks the 
first-order solid-superfluid boundary, the full lines denote the 
continuous transitions between the superfiuid and the fully 
polarised state. 


heavy 6 -electrons and light a-electrons Vab = —Vba 7 ^ 0 . 
Using self-consistent mean-field (BCS-like) theory they 
found solutions with spontaneous on-site hybridisation - 
excitonic condensate. They showed that in case that the 
a and 6 orbitals are of opposite parity the excitonic con¬ 
densation gives rise to a ferroelectric polarisation. Im¬ 
portantly, their excitonic condensate existed also in the 
limiting case of Vab = Vba = 0 suggesting a new ground 
state of the original FKM. This result was challenged by 
subsequent theoretical studies^°“^^ of FKM in one and 
infinite dimension, which found no spontaneous hybridi¬ 
sation. It was pointed out that the method of Porten- 
gen and collaborators missed competing ordered states 
and the possibility of phase separation. Freericks and 
Zlatic^'^ argued that spontaneous hybridisation in FKM 
breaks the local U{\) gauge symmetry, associated with 
the phase of the heavy electron, and thus is prohibited 
by Elitzur’s theorem^®. The question is, however, quite 
subtle as can be illustrated, for example, by the exci¬ 
tonic susceptibility in d = 00 FKM, which has a loga¬ 
rithmic singularity at T = 0^®. The insight into the issue 
can be provided by going into EFKM. The mapping be¬ 
tween 44 > 0 and tatb < 0 models (on bipartite lattice) 
implies that if there is an (ferro-) ordered state in the 
limit tb —7 0 “ there must an antiferro- order in the limit 
4 —7 O'*" and thus FKM is an unstable point between 
these two phases. 


A. Falicov-Kimball model 

The interest in excitonic condensation in FKM started 
with the work of Portengen et a6^®d9 studied FKM 
model augmented with the hybridisation between the 


B. Extended Falicov-Kimball model 

1. Strong-coupling limit 

An alternative way to EC order in FKM is to make 
both fermion species mobile, i.e., introducing EFKM. 



4 


Batista^ studied EFKM at half-filling in the strong cou¬ 
pling limit U ^ ta,tb as an alternative way to electronic 
ferroelectricity. Introducing pseudo spin variables 


t!" = alb^+bla,, rf = iibja, - ajb,), 

< = 4 ®* - bjb,, 


( 2 ) 


he arrived at a strong-coupling effective Hamiltonian, 
which in the case Vab = Vba = 0 reads 


T-LeS = T 


^ II (^11 1 

(a) 


vE’-' 


( 3 ) 

"Heff acts in the low-energy space built of states with 
singly occupied sites only. The coupling constants to the 
lowest order in t/17 read and Kj_ = La¬ 

belling the two local states | fi) = oJlO) and | 1^) = 5||0) 
the structure of the S' = 1/2 XXZ-model in an exter¬ 
nal field along the z-axis is apparent. Note that while 
K± can be both positive and negative, is always pos¬ 
itive. There is a well-known exact mapping of (3) onto 
the model of spinless hard-core bosons^® 


iies = e ^ rii + K± ^ (^4 dj + H.c.^ -I- 2X|| ^ mrij , 

* (b> (b) 

(4) 

with e = A— (z is the number of nearest neighbours). 
The operator 4 (d ^) creates (annihilates) a boson on site 
i, and rii = 4dj is the local density operator. The phys¬ 
ical states are constrained to those containing zero or 
one bosons per site. In this language, the | 1^) is the 
local bosonic vacuum and | fi) = 41 '1'*) i® ^ state with 
one boson. The transverse spin coupling translates into 
bosonic hopping, the longitudinal coupling into nearest- 
neighbour repulsion and the magnetic field into bosonic 
chemical potential. Moving between the spin and boson 
formulation has been traditionally used to allow conve¬ 
nient treatment of various models^®’®®. 

The model (3,4) is much studied in the context of cold 
atoms on optical lattices. Its T = 0 phase diagram for 
a square lattice is shown in Fig. 2. Besides the trivial 
phases obtained for large |A|, which correspond to or¬ 
bitals of one flavour being filled and the other empty 
(saturated spin polarisation along the z-axis), there are 
two more phases. The solid phase (z-axis Neel antifer- 
romagnet) favoured by Ary, and the superfluid phase {xy 
magnetic order) favoured by K±. The solid phase is char¬ 
acterised by checker-board arrangement of sites with oc¬ 
cupied a and b orbitals. The solid phase is connected 
to the ground state of half-filled FKM®^. The superfluid 
phase is characterised by a finite expectation value (di), 
i.e., can be described as a condensate of the d-bosons. In 
the language of the EFKM this means a finite {a\bi) ex¬ 
pectation value characterising the excitonic condensate. 
To understand the phase diagram in Fig. 2 one may start 
form the familiar point of Heisenberg antiferromagnet, 
A = 0, K±^ = K\]^. Upon application of a magnetic field 


(—A) the ordered moments pick an arbitrary perpendic¬ 
ular orientation with a small tilt in the field direction 
(superfluid). For \K±\ > the in-plane order is ob¬ 
tained without the external field, while for iK'j.l < itTy a 
finite field is needed to destroy the Ising (solid) order. 

The solid and superfluid phases have quite different 
properties. The solid phase breaks the discrete transla¬ 
tional symmetry. The superfluid phase breaks continuous 
U{1) symmetry associated with the phase of d-boson. For 
K± > 0 the superfluid also breaks the translational sym¬ 
metry with the phase of {di) varying between sublattices. 
Therefore in two dimensions only solid long-range order 
exists at finite temperature®^, while the superfluid phase 
has the form of Kosterlitz-Thouless phase®®’®^. The mis¬ 
match between the symmetries of solid and superfluid 
phases implies a first-oder transition between them or ex¬ 
istence of an intermediate supersolid phase where both 
the orders are present. The stability of the supersolid 
phase is a much studied question in the context of model 
(4) and its generalisations. Investigations on cubic lat¬ 
tice in two^^ and three®® dimensions found the first-order 
transition scenario to be realised. However, a robust su¬ 
persolid phase was found on the triangular lattice®®^®®. 
Existence of a supersolid on triangular lattice is related 
to the fact that the solid and the superfluid adapt differ¬ 
ently to the geometrical frustration. For further reading 
on hard-core boson we refer the reader to specialised lit¬ 
erature. 

We conclude this section by considering the effect of 
cross-hopping. Non-zero cross-hopping Vab, Vba in (1) 
breaks the C/(I) symmetry of EFKM and generates ad¬ 
ditional terms in Hamiltonian (4)®. These are of two 
types. First, the ’correlated-hopping’ terms of the form 
(dj -I- dl){2nj - 1) and of the order {Vabta - Vbatb)/U- 
Summed over the nn sites their contribution can be split 
into a source field for d-bosons, proportional to the mean 
boson density (n), and a coupling of dj -I- d| to the den¬ 
sity fluctuations. With a finite source field the situation 
is analogous to a ferromagnetic transition in an external 
magnetic field, i.e., the sharp transition is smeared out 
and diminishes completely if the external field is compa¬ 
rable or larger than the internal (Weiss) fields. In some 
symmetries the contributions of individual neighbours to 
the source fields add up to zero, e.g., in the case of a- and 
&-orbitals in (1) of different parity®. In this case, phase 
transition is still possible. However, the second type of 
terms generated by cross-hopping of the form d^dj and 

d|^d] and of the order VabVba/U reduce the symmetry of 
the system to Z 2 . Depending on the sign of this term 
the order parameter (d^ = {a\bi) is either real or imag¬ 
inary. The meaning of the phase of the order parameter 
{a\bi) was discussed by Halperin and Rice®®. In case 
or real orbitals ■i/'o(r) and "(pbi^), real (a-&i) gives rise to 
charge density modulation while imaginary (a|6i) gives 
rise to periodic current pattern. In a model with inver¬ 
sion symmetry about the lattice sites, studied in Ref. 5, 
Vab = —Vba leads to real (albi) implying an Ising-like 
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ferroelectric transition. Solution with purely imaginary 
(a|6i) was reported in d = 1 EFKM by Sarasua and 
Continentino^® for somewhat artihcial choice of purely 
imaginary Kb = Ka- 

2. Intermediate and weak coupling 

Studies using various techniques and in different di- 
menssionalities, reviewed below, lead to the conclusion 
that the behaviour observed in the strong-coupling limit 
of EFKM extends to the intermediate coupling and con¬ 
nects to the weak-coupling regime. To be able to compare 
these results we have to understand how they are inter¬ 
preted, in particular in d = 1, 2. Following the Mermin- 
Wagner theorem^^, the superfluid, which breaks contin¬ 
uous symmetry, is characterised by long-range order at 
finite T in d = 3, by long-range order at T = 0 and 
algebraic correlations at finite T in d = 2, and by al¬ 
gebraic correlations at T = 0 in d = 1. In the Ising- 
like solid phase a long-range order exists at finite T in 
d = 2,3, while in d = 1 there is a long-range order at 
T = 0. The onset of algebraic correlations was used as a 
criterion to define the transition to the superfluid phase 
in the d = 1,T = 0 Monte-Carlo studies^^. Mean-field 
techniques on the other hand do not distinguish dimens- 
sionalities and long-range order is obtained even if it does 
not exist in the exact solution. In these cases one com¬ 
pares the onset of algebraic correlations in more rigorous 
methods with the onset of mean-field long-range order. 

Using the constrained path Monte-Carlo method 
Batista et obtained the T = 0 phase diagram of the 
d = 1, 2 model shown in Fig. 3. In fact, Farkasovsky^^ 
showed that self-consistent Hartree-Fock method repro¬ 
duces the Monte-Carlo phase diagram in d = 2 remark¬ 
ably well (see Fig. 3) and extended the study to d = 3 
and 4^1. A similar HF phase diagram was ob¬ 
tained by Schneider and Czycholl"^^ for semielliptic den¬ 
sities of states (d = 00 Bethe lattice). In both studies, 
the excitonic phase persists in the limit \tb\ —t 0 for cer¬ 
tain range of A. Due to the symmetry of EFKM under 
tb o —tb, tb — 0 is an unstable fixed point between the 
ferro and antiferro excitonic phases. Interestingly, for 
tb ^ 0 Farkasovsky'*^ finds also a phase where both the 
solid and excitonic orders are present simultaneously^"' - 
a supersolid phase. This result has not been conhrmed 
by other studies yet and the supersolid phase was shown 
to be unstable in the strong-coupling limit^^. Moreover, 
existence of a phase with finite {a^b) at tb = 0, which 
smoothly connects to both tb > 0 and tb < 0 violates the 
Elitzur’s theorem^®. 


3. BCS-BEC crossover 

The similarity of the strong- and weak-coupling phase 
diagrams opens an interesting question of how the physics 
described by the the hard-core bosons (3) evolves into the 



FIG. 3. The T = 0 phase diagram of d = 2 EFKM for 
U/ta = 2 obtained with constrained path Monte-Carlo""^ 
(top) and Hartree-Fock approximation"'^ (bottom). The a 
and e (violet shaded) phases correspond to filled b and a band, 
respectively. The solid 7 phase (grey shaded) is surrounded by 
the excitonic (superfluid) phase /3. The supersolid j3' was not 
found in the Monte-Carlo simulations on square lattice. The 
dotted lines in the top panel represent the strong-coupling 
phase boundaries of Fig. 2. The bottom panel was adapted 
with permission of author from Ref. 42. Copyrighted by the 
American Physical Society. 


Hartree-Fock physics of fermions in (1). Similar question 
is known as the BCS-BEC crossover in the context of su¬ 
perconductivity or crossover between Slater and Heisen¬ 
berg antiferromagnet. The issue of BCS-BEC crossover 
in the excitonic phase of half-filled EFKM was stud¬ 
ied by several authors using random phase approxima¬ 
tion (RPA)"'®, slave-boson'^®’"'^, projective renormalisa- 
tion"'®, variational cluster"'®, exact diagonalisation®®, and 
density-matrix renormalisation group®' techniques. The 
different methods provide quite consistent picture of the 
phase diagram with the provision that long-wavelength 
fluctuations of the order parameters are ignored by some 
of the mean-field methods, which therefore describe crit¬ 
ical phase in d = 1,T = 0 and the Kosterlitz-Thouless 
phase in d = 2 as phases with true long-range order. 

The BCS vs BEC question can be formulated as: Are 
there long-lived excitons above K that condense at the 
transition? The typical T — U phase diagram of half- 
filled EFKM with exciton condensate is shown in Fig. 4. 
Following Ref. 45 we introduce a boson creation opera¬ 
tor d)j = ak_|_q&k and the corresponding propaga¬ 

tor Xq®(a;) = ((dq; d(j))c^. The excitonic transition pro¬ 
ceeds differently on the semi-metal and semi-conductor 
sides. In semimetal, the particle-hole continuum extends 
to zero energy in a finite q-region of the Brillouin zone. 
While long-lived excitons may exist in some other parts 
of the Brillouin zone (Fig. 5b), they are not the lowest en¬ 
ergy excitations and the excitonic transition follows the 
BCS scenario. In semiconductor, the excitonic band in 
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FIG. 4. The phase diagram of d = 2 half-filled EFKM for 
A = 2.4, ta = 1, ff, = —0.8 obtained with RPA (solid lines) 
and slave bosons (dashed lines). Note that the temperature 
is scaled to the maximum critical temperature Tmax, where 
Tmax = 0.361 (RPA) and Tmax = 0.256 (slave bosons). The 
points marked 1-4 are discussed in the text. Adapted with 
permission of authors from Ref. 45. Copyrighted by the Amer¬ 
ican Physical Society. 


the two-particle spectrum extends throughout the whole 
Brillouin zone (Fig. 5e-g) and lies below the particle-hole 
continuum. Bqiuj) , the coherent part of ImXq^(<j-'), can 
be approximated as 

^q(w) = -7rZx(q)(5(a; - a;x(q)), (5) 

where Zx(q) is the weight of the excitonic quasiparticle 
and a;x(q) is its dispersion. Note that both Zxiq) and 
a;x(q) depend on temperature. The number of excitons 
with crystal-momentum q is given by 

-^x(q) = Zx{q)nBE{uJx{ci)), (6) 

where nBE(w) is the Bose-Einstein (BE) function. The 
excitonic transition in a semiconductor is connected with 
the minimum of a;x(q) going to zero at Tc and the re¬ 
sulting divergence of iVx(q), Fig. 6. 

The semimetal and semiconductor side of the phase di¬ 
agram appear to be qualitatively different. It is therefore 
instructive to see that the evolution from a semimetal to 
a semiconductor with increasing U is in fact smooth al¬ 
though they differ by the presence of the excitonic band. 
The key observation is that the parts of the exciton band 
that split off the particle-hole continuum carry little spec¬ 
tral weight Zx (q). Moreover, the notion of infinitely 
sharp Bqiuj) is an idealisation, which holds only for ex¬ 
citonic peak well separated from the particle-hole con¬ 
tinuum. Fig. 5 shows the evolution of Zx{q)- Start¬ 
ing from a semimetal with no exciton band, through a 
semimetal with exciton band in a part of the Brillouin 
zone, the region of ^x(q) = 0 shrinks to a point at the 
semimetal/semiconductor crossover. On the semiconduc¬ 
tor side, Zx(q) remains finite over the entire Brillouin 
zone. Vanishing of the excitonic insulator for large U is 



(a)t/ = 0.50, 
T/Tmax = 1.108 



Z(a)x. q) 

1 

0.6 

0 


(b)f/ = 3.08, ic)U = 5.07, 
T/Pmax = 0.679 T/Pnax = 1.108 





Z(a)x. q) 



{g)U = 50.0, 

T/Pnax = 1.108 


FIG. 5. Evolution of the exciton quasiparticle weight ^x(q), 
obtained with RPA, across the semimetal-semiconductor 
crossover. On the metallic side (a-c) non-zero Zx(q) exists 
only in a part of the Brillouin zone. On the semiconductor 
side (e-g), Zx(q) is finite everywhere in the Brillouin zone and 
approaches Zx(q) = 1 in the strong-coupling limit (g). Panel 
(a) corresponds to a weakly coupled semimetal with no exci¬ 
tonic quasi-particle residues anywhere in the Brillouin zone. 
Panels (b)-(f) correspond to the points 1-4 of Fig. 4. Adapted 
with permission of authors from Ref. 45. Copyrighted by the 
American Physical Society. 


Nx(q) 


■3, Nx(q) ■ 
10^ 10“® 






FIG. 6. The occupancy of the excitonic states Nx(q) at 
points 3 (left panel) and 4 (right panel) of Fig. 4 on the 
semiconductor side of the phase diagram. The diverging oc¬ 
cupancy of the lowest mode, q = (0,0), at P, is clearly visi¬ 
ble. Adapted with permission of authors from Ref. 45. Copy¬ 
righted by the American Physical Society. 


the consequence of uix (q) > 0 being finite for all q down 
to r = 0. 

Another question worth understanding is the relation¬ 
ship of the present picture to the strong-coupling limit 
ta,b ^ U. We point out that to reach this limit one can¬ 
not keep A fixed but have to scale it as 1/U (see also 
Fig. 2). Although the bosonic excitons exist in a semi¬ 
conductor away from the strong-coupling limit, there is 
no separation of the bosonic and fermionic energy scales. 
This is reflected in the Zx (q) substantially different from 
1. The dynamics of the excitons therefore cannot be 
described by purely bosonic Hamiltonian. One has the 
option to keep the fermions in the model, use a more gen¬ 
eral action description of the exciton dynamics or resort 
to an effective bosonic Hamiltonian describing only the 
vicinity of the minimum of the exciton band. While the 
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particle-hole gap grows linearly with increasing U, the 
exciton band remains located around A with its width 
approaching the tah/U scaling. If the particle-hole gap 
is sufficiently large, the excitons are no more dressed with 
the particle-hole excitations and thus Zx{o^ ~ 1 every¬ 
where in the Brillouin zone. The exciton dynamics in this 
limit is described by purely bosonic Hamiltonian (4). 


III. TWO-BAND HUBBARD MODEL (S=l/2 
FERMIONS) 

The two-band Hubbard model (2BHM) generalises 
EFKM (1) to include electron spin. In absence of ex¬ 
ternal magnetic field and spin-orbit coupling the one- 
particle part Ht consists of two identical copies of the 
corresponding terms in EFKM. The interaction part i7int 
is richer. 

772BH =Ht -f Hint 

+ H.c. 

<b>,o- 

+ Y1 + H.C. 

<b>,o- 

+ (7) 

i,(T ^ ' 

Hint =u ^ + U' ^ 

i 

-J {Kan\a + 4^a^sb\sb^a) 

IG 

+H-C., 

i 

with the notation a = —a. There are two general set¬ 
ups in which EC in 2BHM have been studied: a lattice of 
two-orbital atoms and a bi-layer Hubbard model system. 

In the bi-layer model the orbitals at and bi are as¬ 
sumed to be spatially well separated and the exchange- 
interaction is vanishingly small J', J ~ 0. This is equiv¬ 
alent to the dominant-term approximation used exten¬ 
sively for long-range interaction^^. Similarly the on-site 
and inter-site inter-layer tunnelling, Vab, Vba ~ 0, is neg¬ 
ligible due to the spatial separation of layers. 

In two-orbital atoms, the overlap of (orthogonal) a and 
b orbitals gives rise to a sizeable ferromagnetic Hund’s 
exchange J and pair hopping J'. The on-site a-b hop¬ 
ping vanishes either by symmetry or can be eliminated 
by a basis transformation. The cross-hopping 14b, Vba 
is in general non-zero, however, in materials with high 
symmetry it may vanish as well.®^. 


A. Normal state 

Before discussing the ordered phases we briefly sum¬ 
marise the basic properties of half-filled 2BHM without 


-1-1— 
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J=U/4 
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|b> . % . .f. 
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FIG. 7. Left: The cartoon representation of the atomic low- 
spin (LS) and high-spin (HS) states. Right: The phase dia¬ 
gram of d = oo 2BHM with symmetric ta = tb bands (semi- 
elliptic density of states) for J' = J = 11/4-, U' = U — 2J 
(red) and J' = J = Q\ U' = U (blue) obtained by Werner and 
Millis®®. The dotted red line marks the HS-LS crossover. The 
vertical blue line for J — 0 corresponds to a model with 6-fold 
degenerate atomic ground state. U and A are expressed in 
the units of bandwidth W. 


broken symmetry. Its physics at strong and intermediate 
coupling is controlled by competition between the Hund’s 
coupling J and the crystal-field splitting Large 

A favours the singlet low-spin (LS) state, while large 
J favours the triplet high-spin (HS) state, see Fig. 7. 
For finite J the phase diagram, shown in Fig. 7, con¬ 
tains three regions: HS Mott insulator connected to the 
limit U ^ W, with W standing for the bare bandwidth, 
and J ^ A, LS band insulator connected to the limit 
A ^ W,U and a metal connected to the non-interacting 
limit and A < W. The low-energy physics deep in the 
Mott phase is described by S' = 1 Heisenberg model with 
anti-ferromagnetic interaction. The band insulator far 
away from the phase boundaries is a global singlet sepa¬ 
rated by large gap from the excited states. In the vicinity 
of the HS-LS crossover both LS and HS states have to be 
taken into account. The physics arising in this parameter 
region is the subject of subsequent sections. 

The physics of the J, J' =0 model is different. This 
setting corresponds to a bi-layer system where one would 
typically choose U > U'. In this case, the low-energy 
physics for IT ^ [7 and A ^ (U—U') is described by ^ = 
1/2 bi-layer Heisenberg model, while for A ^ {U — U') 
one obtains a band insulator. The region A « (C/ — U') 
is described by exciton-t-J model, discussed later. The 
choice U = U', shown in Fig. 7, is anomalous in the 
sense that the bi-layer Heisenberg region is absent and 
the region A « 0 corresponds to exciton-t-J model with 
two exciton flavours. 

The low-energy physics of the metallic phase is sen¬ 
sitive to Fermi surface nesting, in particular for weak 
coupling. Nesting plays an important role in some popu¬ 
lar models, e.g., the Fermi surfaces derived from a and b 
bands are perfectly nested on cubic lattice both in d = 2 
and d = 3. 







B. Excitonic order parameter 


C. Strong-coupling limit 


Excitonic condensation is characterised by a sponta¬ 
neous coherence between the a- and 6-electrons in (7), 
i.e., appearance of matrix elements of the form 

K7' = {4aha') and = (aL^kV')- (8) 

These can be in general complicated objects charac¬ 
terised by translational symmetry, internal structure and 
spin symmetry. 

We will consider only ordered phases with single-q 
translational symmetry, ~ <5k',k+q- On the square 
lattice we will encounter only ferro-EC, q = 0, and 
antiferro-EC, q = (tt, tt) states, which correspond to uni¬ 
form Fa and staggered Fa ~ (—1)*, respectively. 

The internal structure describes the behaviour of the 
above matrix elements as a function of the reciprocal 
vector k or the relative position — R^/. The internal 
structure reflects the symmetry of the pairing interac¬ 
tion. Isotropic pairing leads to isotropic Jic.k+q and, in 
particular, to a finite local element F^. Since this is the 
case for local Hubbard interaction (7) we can use the 
local element Fa as an order parameter for the exciton 
condensation. The spatial decay of the a — b coherence 
can be quantified by a correlation length rcoh defined as 


was studied by several authors for EFKM^®“®^ as well 
as 2BHM^®. Since the largest contribution to the pair¬ 
ing interaction in these models is on-site the correlation 
length reflects the relative size of the pairing field to the 
bandwidth. Weak pairing yields sizeable J^kk only in the 
vicinity of the Fermi level and one ends with rcoh over 
many unit cells (BCS limit). Strong pairing leads to al¬ 
most constant J^kk and Ccoh limited to a few sites (BEC 
limit). 

Finally, we discuss the spin structure of the EC or¬ 
der parameter, which we divide into the spin-singlet and 
spin-tripels parts 

= I + 4 ■ T.v), (10) 


At half-filling and large U,U' ta,tb,Vab,Vba, J, J' 
the charge fluctuations are strongly suppressed. Simi¬ 
lar to the strong-coupling treatment of EFKM, the low- 
energy physics can be described by an effective model 
built on states containing only doubly occupied sites. 
The virtual excitations to states with singly and triply oc¬ 
cupied sites provide inter-site couplings in the low-energy 
model. Using Schrieffer-Wollf transformation to the sec¬ 
ond order one arrives at expressions of the type /U, 
known from the analogous transformation from single¬ 
band Hubbard to Heisenberg model. By construction, 
the low-energy model does not capture one-particle or 
charge excitations, which take place only at energies of 
the order U. The general strong-coupling expressions for 
the case of symmetric bands ta = h were derived by 
Balents®®. In the following we discuss several special pa¬ 
rameter choices of the strong-coupling model, which are 
studied in the literature. 


1. Large J - general formulation 

For a sufficiently large Hund’s coupling J, the 
low-energy Hilbert space can be constructed from 
the atomic HS and LS states: |1) = a|6||vac), 

|0) = -I-a](6j)|vac), | - 1) = a|6||vac), and 

|0) = 6|6j|vac), where |vac) is the fermionic vacuum. 
The local Hilbert space is further reduced if we assume 
easy axis anisotropy and can drop the |0) state. In 
this case only the density-density part of the Hund’s 
interaction contributes.^^ In the following, we will derive 
the strong-coupling model for this case to demonstrate 
the principle. The generalisations are straightforward 
and can be found in the literature®®’^®. Similar to (2), 
the low-energy Hamiltonian can be formulated in terms 
of pseudospin variables. Following Ref. 56, we introduce 
on-site standard-basis operators®® y™-" with matrix 
elements in the local basis 

( 12 ) 

The effective Hamiltonian reads 




i.e., singlet and triplet cj)^ components are defined by 

"" . (11) 

gg ' 

where r are the Pauli matrices and * denotes the complex 
conjugation. Tensorial character of {^®, and the fact 
that their elements are complex numbers allow numerous 
distinct phases as will be discussed in what follows. Some 
of the phases lead to a uniform spin polarisation while 
others have no ordered spin moments. 


=e E E ^ j) 

i,s {ij),s 


+A'. y {l 


+ T^Tf 


(13) 


where s = ±1 and s = —s. Typical hopping processes 
contributing to are shown in Fig. 8. The pro¬ 

cess (i) lowers the energy of HS-LS pair on nn bond 
relative to the LS-LS pair. Therefore it lowers the 
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(i) 


(ii) 


(iii) 


|a> ..*.4.. 

ib >n t n t 


FIG. 8. Typical nn hopping processes that give rise to cou¬ 
plings in the effective Hamiltonian (13, 15, 17). 


energy e of a single HS site on otherwise LS lattice 
relative to the single atom value of Ehs — Els, and 
contributes to a nn repulsion between HS states 
A similar process between two HS states with oppo¬ 
site s gives rise to the exchange term Kq. The pro¬ 
cess (ii) exchanges the HS and LS states in a nn bond 
and introduces quantum fluctuations in the model. Ne¬ 
glecting the cross-hopping contributions (see Ref. 58 for 
the general expressions including cross-hopping correc- 

tions) the coupling constants read: e = A — 3 J — z , 

rr. — 2tati, _ (f2 I f2\ 2U-U'+2J j _ *o+*b 

J^± — J// , ^||—+ [//([/+J) J ana Ao — [/+j , 

where z is the number of nearest neighbours. Fi¬ 
nally, the process (iii) converts HS-HS pairs with 
zero total moment into LS-LS pairs and vice versa. 
This process is possible only with finite cross-hopping. 


Hamiltonian (13) can be formulated in terms of hard¬ 
core bosons. In this picture, |0) is identified with bosonic 
vacuum and |s) with a state containing one boson of 
flavour s. However, practical implementation of the hard¬ 
core constraint prohibiting more than one boson per site 
is complicated. The standard way to treat the constraint 
is to introduce a new vacuum state |0) and Schwinger- 
like bosons: |s) = |0) = /i'^|H). The physical states 

are required to obey the local constraint 




(14) 


Rewriting (13) in terms of d- and /i-boson, using the re¬ 
placement T*® = d\ and introducing the d-number 
operator rii = d\ ^ T/^ and the spin opera¬ 
tor S'f = X;* one arrives at 

* (b),s 

(15) 

(b) (b) 

-^1 E +id.c.). 


advantages, the introduction of d-boson for |0) treats the 
LS and HS states on equal footing and thus is well suited 
for e < 0, where |0) is not the atomic ground state and 
thus cannot be viewed as the vacuum state. 

The effective Hamiltonian for the case with SU (2) spin 
symmetry has a similar structure, but differs by the pres¬ 
ence of a third bosonic flavour |0) = dj|fl). Introducing 
a cartesian vector with components 

^x\ / ^-1 ~ 4 \ 

4 = ^ *(4 1 +dl) , (16) 

4) V2dl ) 


the effective Hamiltonian can be written in a compact 
form 


= £ E E • 444 + H.c) 

* (b> 

+ E "^*^3 +Ko'^Si- Sj 
{ij) (ij) 

+ iLiE (4 -444 + 

(b> 

(b) 


(17) 


Here, are spin S' = 1 operators with Sf = ^(T/® -|- 
7;°-i+T4+71-10), Sf = ^(Ii0i+7;-i0-Iii0-7;0-i). 
Expressed in vector notation, takes the form of a cross 
product 

Si = —id| A di. (18) 

Unlike the Ising case (15) where cross-hopping is neces¬ 
sary to generate Ki, in the SU (2) symmetric case (17) Ki 
appears also without cross-hopping if the pair-hopping 
term is finite, J' ^ 0. The last term with coupling con¬ 
stant K 2 ~ Vbata + Vabtb, which couples the d-operators 
to the spin operators, does not have an analogy in the 
Ising case. The full expressions for the coupling constants 
can be found in Ref. 58. 

Exciton condensation in (17) is characterised by a finite 
expectation value (d^/ii), related to the spin-triplet order 
parameter (11) of 2BHM by 

(did.) = cl>l/V2. (19) 

Models (15) and (17) for special choices of the coupling 
constants are known under their own names. We discuss 
these cases below. 


2. Blume-Emmery-Griffiths model 


This allows us to interpret the Kl term as nn hopping of 
d-bosons, K\\ as nn repulsion between the d-bosons, Kq as 
nn spin-spin interaction and Ki as pairvise creation and 
annihilation of d-bosons on nn sites. Besides technical 


Without cross-hopping and one fermionic species im¬ 
mobile, tb,Vab,Vba = 0, (15) becomes purely classi¬ 
cal, with K±,Ki = 0. Known as the Blume-Emmery- 
Griffiths (BEG) model®, it was originally introduced to 
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FIG. 9. The phase diagram of BEG model for Tfn/lfo = 3.5. 
Besides the normal phase (N), it contains the solid phase (S), 
antiferromagnetic phase (AFM) and ferrimagnetic phase (I). 
The data were taken from Ref. 60. 

describe mixtures of ^He and "^He. In the standard for¬ 
mulation of the BEG model discrete index s is used to 
describe the local state s = ±1 corresponding to | ± 1) 
and s = 0 is assigned to |0). The model is then written 
as 

Hbeg = +Ko^s,Sj. ( 20 ) 

* (b> (ij) 

The BEG model found its use in many areas of statis¬ 
tical physics. Despite its simplicity it has a rich phase 
diagram®^. The phase diagram of the model on bipartite 
lattice does not depend on the sign of Kq as the ferro¬ 
magnetic Kq < 0 and antiferromagnetic Kq > 0 mod¬ 
els can be mapped on each other by the transformation 
Si — >■ (— l)*Si. For Ffii < \Ko\ there are only magnetically 
ordered (AFM) and normal (N) phases. At T = 0 the N 
phase corresponds to empty lattice (rii) = 0, while in the 
AFM phase there is one boson on each site (n^) = 1. A 
typical phase diagram for ATy > jATol is shown in Fig. 9. 
Besides the magnetically ordered (AFM) phase there is 
the solid (S) phase (called antiquadrupolar in Ref. 60) 
characterised at T = 0 by bosons occupying one sublat¬ 
tice, the other being empty, (n^) = ^(1 — (—1)*). Note 
that this phase has a residual spin degeneracy on the oc¬ 
cupied sites. This is a consequence of the exchange inter¬ 
action beyond nn being strictly zero. In a more realistic 
model one expects the occupied sites to order magnet¬ 
ically at a sufficiently low temperature. Goexistence of 
the magnetic and solid order (I phase in Fig. 9 ) is found 
at finite temperature in a narrow range of e separating 
the magnetic and solid phases. For positive e at T = 0, 
the system is in the vacuum (empty lattice). Interest¬ 
ingly, for moderate e > 0 the solid phase is found at 
elevated temperature. This reentrant behaviour of the 
solid phase was found also for finite 4 in Monte-Carlo 
simulations of the spinless bosons^^ as well as in DMFT 
simulations of 2BHM®^. Generalisation of BEG model to 
the case with spin-rotational symmetry (17) is straight¬ 
forward. On the mean-held level it leads to quantitative 
modihcation of the phase diagram. 


3. Bosonic t-J model 

Next, we discuss models (15) and (17) in the param¬ 
eter range where they describe conserved bosons. We 
assume that both a and b electrons are mobile generat¬ 
ing a hnite hopping K± in (15) and (17) of the d-bosons. 
In addition, we require that Ki,K 2 = 0, which is the 
case for 146 , 14a, J' = 0- Hamiltonians (15) and (17) de¬ 
scribe interacting bosons carrying Ising and Heisenberg 
spin S' = 1, respectively. For suitable parameters, the 
system may undergo a BE condensation characterised by 
a complex vector order parameter </>( (19). In case of 
the Ising interaction (15) the vector <pl is conhned to the 
xy plane. The vector character of the order parameter 
adds additional structure to the condensation of spinless 
bosons discussed in Sec. HBl. 

The role of Kq . The SU (2) symmetric Hamiltonian 
(17) describes spinful S' = 1 bosons with inhnite on-site 
repulsion, nn repulsion Ary and nn spin-exchange Kq. 
Hamiltonians of this kind have been studied both the¬ 
oretically and experimentally for cold atoms in optical 
traps. Interestingly, in a genuine system of bosons with 
spin-independent interactions the structure of the exact 
ground state prohibits the spin-exchange from appearing 
in any low-energy effective Hamiltonian®^. As pointed 
out in Ref. 62, it can only arise as a low-energy effective 
description in fermionic system - as in the present case. 

BE condensation in the continuum version of (17) was 
studied in Refs.®®’®"*. The sign of the exchange was 
shown to play a crucial role for the properties of the 
superfluid phase as it determines the residual symme¬ 
try of the BE condensate. Antiferromagnetic exchange 
(ATo > 0) selects the so called polar state characterised by 
{4>l)* A (/)* = 0, while ferromagnetic exchange {Kq < 0) 
selects ferromagnetic state, where \{4>l)* A cfH is max¬ 
imised. Different residual symmetries of these states re¬ 
sult in different low-energy dynamics and qualitatively 
different behaviour of topological defects (vortices)®®’®"*. 
The polar and ferromagnetic phases are found for Ising 
spins (15)®®, although the topological aspects that deter¬ 
mine the low-energy excitations and possible topological 
defects are different from Heisenberg spins (17). 

Mean-field phase diagram. The continuum model may 
be viewed as an effective description of the lattice model 
at low boson concentrations. At higher boson concen¬ 
trations other phases, e.g., these present in the BEG 
phase diagram Fig. 9, exist on a lattice and compete 
with the superfluid. In Fig. 10 we show the phase di¬ 
agram of (17) on a square lattice {z = 4) obtained with 
a mean-field decoupling of the pseudospin variables (13) 
TiTj « {Ti)Tj Ti{Tj) — {Ti){Tj). In absence of con¬ 
densed bosons, this mean-field theory reduces to that of a 
slightly modified BEG model (Heisenberg instead of Ising 
spin S' = 1) with phase boundaries marked with the red 
lines. The grey area marks the superfluid (exciton con¬ 
densate) phase - the different shades of grey correspond 
to different values of Arj_. The transition between the 
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FIG. 10. Mean-field e-T phase diagram®® of the bosonic 
t-J model (17) for the parameters: Kq = 0.041, 7F|| = 0.164, 
and IFx^l-O.lOS, -0.076, -0.023}. Solid lines mark continu¬ 
ous transitions, while dotted lines mark first-order transitions. 
The red lines correspond to the phase boundaries of the S' = 1 
BEG model with its phases described in Fig. 9. The shaded 
regions correspond to the EC phase for the three choices of 
K± (with the brightest colour for largest \K±\). Red lines 
inside the EC phase have no meaning and are shown only to 
put several phase diagrams in one figure. 


normal and the EC phase is continuous, while the tran¬ 
sitions to the other ordered phases are of the first-order 
as dictated by symmetry. The intermediate phases, such 
as supersolid (with S and EC orders) or AFM-supersolid 
with (AFM and EC orders), which would allow contin¬ 
uous transitions, do not exist or are thermodynamically 
unstable. 

At large e the ground state is an empty lattice, a state 
that is connected to normal Bose gas at elevated tem¬ 
peratures. Upon reduction of e the system undergoes a 
continuous transition to the superfluid EC phase with 
determined by |Ar_L| and e. The sign of K± determines 
the periodicity of the order parameter: for K± < 0 the 
system goes to uniform ferro-EC state (not to be con¬ 
fused with ferromagnetic EC state) with </>* = cf>, while 
for K± > 0 the system goes into antiferro-EC state with 
(pl = (—!)*</>, the models with ±|Ai_L| are connected by 
the gauge transformation di —>■ {—lydi. The BEG phase 
boundaries are not affected by K±. 

Another parameter that determines the nature of the 
EC phase is the exchange coupling Kq. The phase di¬ 
agram in Fig. 10 was obtained for antiferromagnetic 
Kq > 0, but some of its features remain unchanged when 
the sign of Kq is flipped. In particular, the BEG phase 
boundaries are unchanged. This is so because the on-site 
and interaction terms in the Hamiltonian are invariant 
under the exchange of spin species on one of the sublat¬ 
tices {dis —>■ dis), which maps ferromagnetic BEG to an¬ 
tiferromagnetic BEG model. This transformation works 
only for Ising spins, but in the mean-field treatment the 
difference between Ising and Heisenberg spins disappears. 
The transformation changes the hopping term in (15,17) 
and thus the argument cannot be used for the EC phase. 
Nevertheless, the mean-field phase boundary between the 


normal and EC phase does not depend on Kq at all, be¬ 
cause it contributes to the free energy in the order On 
the other hand, the first-order phase boundaries depend 
Kq and its sign. 

Ferromagnetic EC state. It is instructive to see where 
exactly the difference between the ±|Aro| models comes 
from. First, we point out that an order parameter of 
the form cpj = 4> implies a uniform magnetisation 

(Sj) ~ icp* A 4> irrespective of its periodicity q. Second, 
let us consider the free energy of the FMEC and po¬ 
lar EC states with the same magnitude \4>\ of the order 
parameter. The on-site, hopping and interaction terms 
contribute the same for the two states. However, the 
exchange energy in the polar EC state is zero while the 
FMEC state has finite exchange energy. Therefore, like 
in the continuum models, Kq > 0 leads to a non-magnetic 
polar EC state, while Kq < 0 leads to FMEC with a fi¬ 
nite uniform magnetisation. For the same lATol, FMEC/S 
boundary for ATq < 0 is shifted in favour of the EC phase 
compared to the polar-EC/S boundary for Kq > Q be¬ 
cause the FMEC energy is lower than the corresponding 
polar EC energy. 

Continuous transition between the FM and FMEC 
phases is allowed by symmetry. While we are not aware 
of an explicit calculation of the FM/FMEC transition for 
Kq < 0, one can gain insight from the T = 0 mean-field 
wave function, which has the product form^® 

1^) =n, + d.d^) |H), (21) 

with the coefficients fulfilling -|- |(isp = 1. The 
magnetisation (18) in terms of the EC order parameter 
(19) cp = \f2dFh is given by 

(d/|S,|vI/) = -*d*Ad=^<A*A0. (22) 

This shows that at a continuous FMEC/FM transition h 
goes to zero, while the magnetisation reaches smoothly 
its saturation value. 

Two flavour models. We are of not aware of specific 
theoretical studies of spinful hard-core bosons (15,17) on 
a lattice. However, models (15,17) can be viewed as spe¬ 
cial cases of more general multi-component boson sys¬ 
tems. Bosonic t-J model describing mixtures of two bo¬ 
son species with hard-core constraint have been inves¬ 
tigated in several studies. This set-up is similar to the 
model (15) with Ising spin and Ki = 0, which conserves 
both bosonic species separately. In the following, we dis¬ 
cuss briefly how and which of the results of these studies 
can be used to describe (15). 

The bosonic t-J model is usually formulated in terms 
of / = 1/2 pseudospin: = d\d^ — dLid_i, = did_i, 

and I~ = d}_^di, and anisotropic XXZ inter-site ex¬ 
change interaction. While the definition of spin operator 
Sz in (15) coincides with that of 1^, the transverse com¬ 
ponents /'*■ and I~ have no spin counterpart in (15). 
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Therefore a great care is required when interpreting the 
results of bosonic t-J model in the language of model 
(15). In particular, the z-axis pseudospin order (often re¬ 
ferred simply as an anti-ferromagnetic order) corresponds 
to a true spin order in (15). However, the xy pseudospin 
order (e.g. a;i/-ferromagnet of Refs.®^’®®) does not cor¬ 
respond to a spin order in (15) or its 5'[/(2)-symmetric 
generalisation (17). 

Another difference between the bosonic t-J model and 
(15) concerns the hard-core constraint. In most studies 
on the bosonic t-J model simultaneous presence of two 
bosons of different flavour on the same site is allowed ei¬ 
ther explicitly (The hard-core constraint applies only to 
particles of the same species.) or implicitly (Double occu¬ 
pancy is allowed in the high-energy model from which the 
t-Jmodel is derived.), giving rise to the xy pseudospin 
exchange. On the other hand, the total hard-core con¬ 
straint and the absence of the xy exchange in (15) are 
viewed as kinematic constraints®® . This calls for care 
when using the phase diagrams of bosonic t-J model since 
existence of some phases, e.g., the counter-superfluid of 
Ref.®®, depends crucially on the xy exchange and so such 
phases are not present in (15). 

In the following we discuss two types of studies on the 
bosonic t-J model, which can be translated to provide 
information about (15). The first concerns the tran¬ 
sition between AFM and SF phases and the possibil¬ 
ity of intermediate AFM-supersolid phase. Boninsegni^® 
and Boninsegni and Prokofev^^ studied the bosonic t-Jz 
model, corresponding to (15) for K\\ = —Kq < 0 using 
Monte-Carlo simulations. Since the solid phase is absent 
for ATii < one may expect AFM-SF transition 

for moderate Kq (For strong Kq the BEG model pre¬ 
dicts first order transition between the Mott AFM and 
vacuum states at low temperatures.) The Monte-Carlo 
calculations find a first-order AFM-SF transition with 
no indication of AFM-supersolid characterised by coex¬ 
istence of the AFM and SF order parameters. 

Numerical calculations with repulsive [Kq = > 

0), motivated by observation of the supersolid in the 
spinless case, have been reported on a triangular lattice 
recentlyAlthough double occupancy by different 
bosonic flavours was allowed and the lattice differs from 
the square lattice, the basic features common to Fig. 10 
were found in the phase diagram for large on-site (inter¬ 
species) repulsion. In particular, the sequence of T = 0 
phases ’empty lattice-SF-S-SF’-AFM’ with decreasing e 
can be expected. 


4- Bi-layer Heisenberg model 

Another special case of the S'[/(2)-symmetric model 
(17) studied in the literature corresponds to bi-layer 


Heisenberg model 

~ Jx ^ ^ ' ^ib “b (/|| ^ ' ^jm h • ^ 

i im 

(23) 

with S' = 1/2 spin operators i being the site index 
and m = a, & a layer index. The model arises as the 
large U limit of the half-filled bi-layer Hubbard model 
(7), which describes two identical layers, indexed by the 
orbital index m, coupled on inter-layer rungs. This sit¬ 
uation corresponds to the parameter set U \ta\ = 
\tb\, {A,U',J'} = 0, J < 0 in (7), where the anti¬ 
ferromagnetic coupling on the rung J± = J may arise 
from inter-layer tunnelling. Integrating out the charge 
fluctuations one arrives at (23) with Jy = 4t^/[/. Model 
(17) is obtained by going to the singlet-triplet basis, 
which diagonalises the rung exchange (local) part of (23). 
Introducing a map 

Sia,b = ^ (±tjs^ ± A (24) 

one arrives at^’^^ 

Hhi-H = • ti — 3s|^Si ^ 

I 

(b> 

+ YEH)"(tlAt,) • (tJ-At^.) (25) 

(A) 

(A) 


where Si and ti,z) are bosonic operators 

and the physical subspace fulfils the local constraint 
-I- 1 ■ • t ■ = 1. Hamiltonian (25) is equivalent to (17) 
with parameters ATy =^2 = 0, K±_ = Kq = Ki = Jy/2 
and e = Jl. The main difference of (23) to the bosonic 
t-J model consist in the presence of non-zero ATi, which 
leads to the number of d-bosons being a non-conserved 
quantity. While the system undergoes transitions^"* to 
phases characterised by non-zero cf>i = the phase 

of (pi is not arbitrary and thus the transitions cannot be 
viewed as Bose-Einstein condensation of spinful bosons. 
In zero magnetic field h = 0 the ordered phases are char¬ 
acterised by real (pi which breaks the SU{2) symmetry 
of (23) and corresponds to a Neel state with ordered mo¬ 
ments (Sia^b) = '^(pi, see Fig. 11. Non-zero magnetic field 
h 0 reduces the symmetry of (23) to U{1). The vector 
(pi gets oriented perpendicular to h and acquires an imag¬ 
inary part perpendicular to the real one. The real part of 
(pi describes the AFM order, while i(p* A (pi describes the 
FM component along h. Breaking of the U{1) symmetry 
can be described as BE condensation of spinless bosons. 
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FIG. 11. r=0 phase diagram of the bi-layer Heisenberg 

obtained by Sommer, Vojta and Becker. Reproduced from 
Ref. 74 with kind permission of The European Physical Jour¬ 
nal (EPJ). 


This is particularly easy to see for Jj_ Jy, because close 
to the transition h « Jy and model (17) reduces to the 
Hamiltonian of spinless bosons for the flavour polarised 
along the field direction. 

BE condensation in quantum magnets has been an ac¬ 
tive area of research and the reader is referred to the 
recent review^^ for further reading and references. In 
general, the basic difference between BE condensation in 
spin systems and excitonic condensation discussed here 
consists in the microscopic origin of the 17(1) symme¬ 
try broken by the condensate. In spin models, it is the 
spin rotation in the xy-plane of systems with uniaxial 
anisotropy and condensation refers to some kind of in¬ 
plane magnetic order. The 17(1) symmetry in excitonic 
systems is more abstract and refers to the arbitrariness 
of the relative phase between a and b orbitals in (7). 


5. Exciton t-J model 

The strong-coupling limit of 2BHM with anti¬ 
ferromagnetic Hand’s exchange was studied by Rade- 
maker and collaboratorsThey started from the 
antiferromagnetic bi-layer Heisenberg model (25) and 
considered doping one layer by holes and the other 
by the same amount of electrons. This is equivalent 
to introducing the crystal-field A in (7) while keep¬ 
ing the total electron concentration at half-filling. The 
anti-ferromagnetic Hand’s coupling J < 0 (correspond¬ 
ing to inter-layer exchange) was assumed to originate 
from small, but finite inter-layer tunnelling. Starting 
from the undoped Heisenberg limit with on-site states 
|lm) = and |00) = s'^|n) = ~ aj6|)|vac), 

the authors of Refs. 76-80 introduce a bound exciton 
state |0) = ft.^|n) to the system. The new Hamiltonian, 


which is a special case of that discussed in Ref. 56, reads 

H = iJbi-H + M XI ^ H 

i (ij) 

/ X ( 26 ) 

+ <ex X + H.c.j , 

(b) 

with the hard-core constraint hlh^ +sjsj -I- tj • t ■ = 1. 
The exciton site energy fj, controls the exciton concentra¬ 
tion p. The spin-singlet bosons h and s appear symmet¬ 
rically in (26), except for the term flffs s which does 
not have a tH^h h counterpart. The exciton hopping 
tex and repulsion V play roles analogous to K± and RTy , 
respectively, in (15, 17). 

Note that the notion of exciton and vacuum is ex¬ 
changed with respect to that introduced in Sec. HI Cl. 
Unlike the unambiguous vacuum state of real bosons, the 
meaning of vacuum for hard-core bosons is ambiguous, 
similar to the notion of electron and hole for fermions. 
The condensation of hard-core boson generally means 
that the system is in a quantum-mechanical superpo¬ 
sition of the vacuum and one-particle states, whatever 
their definition is. 

In Refs. 76 and 77 Rademakeret al. studied the prop¬ 
agation of a single |0) exciton in the model with antifer¬ 
romagnetic coupling J||, Jj^ > 0. The propagation of ex¬ 
citon strongly depends on the ratio of the Hund’s (inter¬ 
layer) and intra-layer exchange J_l/J||. For Jl J\\ 
the ground state of undoped system is a product of lo¬ 
cal singlets |00) and the exciton can propagate as a free 
particle forming a band with the width 2ztex- In the 
opposite limit, J± <C J\\ the system consists of weakly 
coupled AFM layers. It is well know from the fermionic 
t — J model that the motion of a hole in the AFM back¬ 
ground is strongly inhibited since a moving hole disturbs 
the AFM order. This physics is also reflected in the mo¬ 
tion of an exciton. For Ux ^ <7||, exciton propagation 
is severely limited and the exciton bandwidth is reduced 
to the order tex/'^||- In the opposite limit tex >7||, the 
exciton can explore neighbouring sites over larger dis¬ 
tances which gives rise to the typical incoherent string- 

state spectrum^®77,80 

Exciton condensation in the excitonic t-J model was 
studied in the Refs. 78 and 80. The general features of 
the mean-field phase diagram, shown in Fig. 12, are sim¬ 
ilar to the case of antiferromagnetic Hund’s exchange in 
Fig. 10 and can be traced back to their ’common an¬ 
cestor’ in the XXZ model. In particular, there are three 
basic phases: AFM, solid and the EC superfluid and first- 
order AFM/solid, solid/superfluid and AFM/superfluid 
transitions. We can compare these Figs. 12 and 10 keep¬ 
ing in mind that increasing e in Fig. 10 corresponds to 
increasing p in Fig. 12. In both models the T = 0 
solid (S) phase is absent for large exciton hopping 
(tex)- At intermediate hopping K_i_ (Ux), the AFM- 
EC-S-EC sequence of phases is found in both mod¬ 
els, while at small hopping there is a direct transition 
between the AFM and S phases. The small positive 
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FIG. 12. The T — 0 phase diagram of the exciton t-J model 
for J± — 0.02, J|| = 0.125 and V = 2. The parameters t 
and p denote the exciton hopping amplitude and concentra¬ 
tion, respectively. The names of phases are the same as in 
Fig. 10: antiferromagnet (AFM), exciton condensate (EC) 
and solid (S). The slash, e.g., ’EC/S’, denotes phase sepa¬ 
ration. Adapted with permission of authors from Ref. 80. 
Copyrighted by the American Physical Society. 

Jj_ selects the spin-singlet condensate which is charac¬ 
terised by spontaneous coherence between the |00) and 
|0) states. The mean-field ground state can be written 
as 141) = \li{yfph\ + yjl — /9s-)|r2). An interesting conse¬ 
quence of the singlet condensation is an enhanced prop¬ 
agation of spin excitations (triplons). The authors of 
Refs. 78 and 80 observed that the triplet excitations in 
the condensate of mobile excitons tex ^ J\\ propagate 
faster that in the quantum paramagnet. Surprisingly, 
they found the effective triplon hopping scales with the 
density of the condensate psv = -y/p(l — p) rather than 
the exciton density p. 


D. Weak coupling 

In the weak-coupling (BCS) limit the formation of ex¬ 
citons and their condensation take place at the same tem¬ 
perature. The physics of the system in this limit can be 
described by approaches such as Hartree-Fock approxi¬ 
mation and RPA. The key feature of the weak-coupling 
theories is that EC instability is driven by nesting be¬ 
tween the Fermi surface sheets formed by the valence 
and conduction bands. The lack of perfect nesting in 
real materials is likely one of the reasons why EC is rarely 
found in nature. The weak-coupling theory of the exci- 
tonic condensation in systems with equal concentration of 
holes and electrons was developed in 1960’s®’^° and sum¬ 
marised in the review articles of Halperin and Rice^^’®^. 
The pairing glue considered in the weakly-coupled semi¬ 
metals or semiconductors comes from the long-range part 
of the Coulomb interaction. The exchange part of the 
long-range Coulomb interaction is small, i.e. similar to 
the choice J, J' « 0 in (7), and one has to consider both 
spin-singlet and spin-triplet pairing. Halperin and Rice^^ 



FIG. 13. The mean-field phase diagram of doped 2BHM as a 
function of external field h and the chemical potential p (left) 
and as a function of doping x (right). The solid lines in the 
left panel mark first-order transitions. The light blue areas 
(marked D) in the right panel represent phase coexistence 
regions. The meaning of the different phases is described in 
the text. Adapted with permission of authors from Ref. 85. 
Copyrighted by the American Physical Society. 


classified the possible excitonic condensates in systems 
with single Fermi surface sheet per band into the charge- 
density wave (real singlet), charge-current-density wave 
(imaginary singlet), spin-density wave (real triplet), spin- 
current-density wave (imaginary triplet) type®^. In case 
that there are multiple Fermi surface sheets related by 
point-group symmetries, e.g., as in hexaborides, a more 
complex symmetry classification is necessary®®’®"^. The 
classification of Ref. 11 includes only the non-magnetic 
polar EC states. In mid 1970’s Volkov and collaborators 
showed that in a doped material with unequal number 
of electrons and holes ferromagnetic EC state develops 
characterised by simultaneous presence of finite singlet 
and triplet components 1^-14 _ 

We briefly review the mean-field (Hartree-Fock) theory 
of (7) and consider the simplest case of uniform EC order 
{tatb < 0). The mean-field Hamiltonian that allows polar 
as well as ferromagnetic EC order reads 

77 — ^ ^ (eka*^CT{7^ h(j • ^kfT^k(T^ 

k,cr(7' 

+ i^kbSaa' — ii-b ■ Tcra') bl^cr\a' (27) 

k,<7(7' 

~ 'y ^ ^(As(5o.(t' + At • Xo-o-') -\- H.c.^ . 

k,<7(7' 

Here, ako-’^ko- Fourier transforms of respec¬ 

tively. The crystal-field splitting A as well as the spin- 
independent part of the self-energy are absorbed in band 
dispersions eka,ekb- In the following discussion we will 
assume a to form the conduction band and b to form the 
valence band. The Weiss fields ha,ht,As, At are given 
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by 

ha = h + ^ (aLaka')'^^<^' 

kjCrcr' 
jjf _^ 

~ 2]v ^ (^kcT^kcr'}^'^'^' (28) 

k.,(7<7' 

Jjf _^ 

~ 2iV ^ (®k(T^kCT')'’’aa', 
k.crtr' 

where h is the external magnetic field (acting on spin 
only). The field h;, is defined as ha with the orbital 
flavour replaced. Note that for model (7) with local in¬ 
teractions {As, At) = where (j)s, (pt are the 

local order parameters. The generalisation to models 
with non-local interaction, which leads to k-dependent 
Weiss fields, is straightforward and can be found in the 
literature. 

Bascones et al.^^ studied model (27) as a function of 
doping and external field h at T = 0. The phase dia¬ 
gram, shown in Fig. 13, contains four phases: the normal 
phase (N), polar excitonic insulator (El) phase and two 
metallic FMEC phases called NC and COL in Ref. 85. 
To illuminate the nature of these phases it is helpful to 
introduce A^a' = ^{As6aa' + At ■ Taa')- 

In the zero-field El phase, the singlet and triplet EC 
orders are degenerate. A finite field h, assuming h = 
hzi and > 0 in the following, lifts the degeneracy. 
The undoped system selects a triplet state with A-|-|- = 
= 0 and |A.|.^| > |A^.|.|. For sufficiently large the 
system enters the EI2 phase of Fig. 13 with A^^ = 0. 
Upon doping two distinct FMEC phases are found. In 
the NC the only non-zero element of the oder parameter 
is A.^^ 7^ 0. The EI2 and NC phases are distinguished by 
presence of a h-dependent gap between the uncondensed 
and 67 bands. The order parameter in the NC phase 
is purely spin-triplet and the phase coincides with the 
EMEC phase of S' = I bosons in the strong-coupling 
limit. In the COL phase, the only non-zero element of 
Aacr' is A^^ 7^ 0. In this phase, predicted by Volkov et 
al}^ the singlet and triplet order parameters mix with 
equal weight (A^ = —Af). 

The h = 0 spin-wave spectrum of the COL state for a 
general chemical potential fi is characterised by two gap¬ 
less modes with a quadratic dispersion at small |q| and 
an additional soft but gapped mode®^. At a special value 
of /X the soft mode becomes gapless and the spectrum 
has one quadratic and two linear modes. The mean-field 
phase diagram and magnetic excitations of the polar El 
phase in the undoped model were studied by Brydon and 
Timm®® and Zocher et at.^"^ using RPA. They observed 
acoustic-like modes with linear dispersion at small |q| 
predicted by Kozlov and Maksimov®® and Jerome et 
Moving away from the strong-coupling limit the sepa¬ 
ration of energy scales of exciton formation and exciton 
condensation is progressively less well defined and even¬ 
tually these scales are not separated at all. As in the spin¬ 
less case of EFKM, EC exists also at weak-coupling and 


one can follow the BEC-BCS crossover as the interaction 
strength is lowered. The weak-coupling methods, such 
as Hartree-Fock approximation and RPA, were applied 
to study the EC in its early days and are summarised in 
the review article of Halperin and Rice^^’®^. 


E. Intermediate coupling 

Investigations of systems with intermediate-coupling 
strength are notoriously difficult due to the lack of small 
parameters. The general approaches to this problem 
include numerical simulations of finite systems such as 
exact diagonalisation or quantum Monte-Carlo (QMC) 
methods, large-N expansions, and embedded impurity or 
imbedded cluster methods such as dynamical mean-field 
theory®® (DMET), variational cluster approximation®^, 
dynamical cluster approximation or cluster DMET®®. A 
major obstacle in simulation of ordering phenomena with 
finite-system methods is the necessity of scaling analy¬ 
sis, i.e., the cluster size must be large enough to show 
the ’diverging’ correlation length. Monte-Carlo simula¬ 
tions on large clusters are available for many bosonic and 
spin systems, but usually not for fermions. Rademaker 
et af^ applied the determinant Monte-Carlo method to 
bi-layer Hubbard model and were able to demonstrate 
an enhanced response to the excitonic pairing field, but 
could not reach temperatures below Tc- Besides the 
Green’s function QMC methods for calculation of cor¬ 
relation functions, wave function QMC approaches can 
be applied to variational search for ground states. While 
we are not aware of variational QMC studies of excitonic 
condensation in Hubbard-type lattice models, variational 
QMC has been used to study the corresponding contin¬ 
uum problem®®. 

The DMET methods have been very valuable for in¬ 
vestigation of Hubbard model and its multi-band gener¬ 
alisations in the past two decades. However, most appli¬ 
cations of DMFT so far focused on one-particle quanti¬ 
ties and normal (paramagnetic) phase. With an excep¬ 
tion of a multi-band study of Ref. 94, linear response 
calculations needed to identify instabilities have been so 
far limited to simple one-®® and two-band models®®’®®’®^. 
Recent DMFT investigations of ordered phases in mod¬ 
els as simple as 2BHM uncovered a rich physics. Capone 
et a/.®®’®®, Koga and Werner^®® and Vanhala at al.^®^ 
found various forms of superconductivity in 2BHM at 
half-filling. Recently, dynamical cluster approximation 
was also applied to study excitonic condensation.®®® 


1. Half filling 

Kunes and collaborators studied 2BHM (7) in the 
vicinity of spin-state transition, see also Fig. 7, using 
DMFT®®’®®’®®>®®2 for 2BHM (7) with Vab,Vba,J' = 0 
and density-density interaction. In Ref.®®, they reported 
observation of the solid phase for the case of strongly 



16 




FIG. 14. Left: The instability temperature of the normal 
phase for various band asymmetries tbjta and crystal-fields A. 
The squares mark the divergence of the orbital susceptibility 
(instability towards the solid order), the circles correspond 
to the divergence of the excitonic susceptibility. The sum 
ta + tb kept fixed. Data taken from Ref. 58. Right: The 
magnitude of the order parameter \(p{T)\ along the trajectory 
marked by arrow in the left panel. The specific heat per site 
along the same line is shown in the Inset. Adapted from 
Ref. 102. Copyrighted by the American Physical Society. 


asymmetric bands \ta\ 3> |t&|. The normal/solid phase 
boundary qualitatively agrees with strong-coupling BEG 
model. In particular, the re-entrant transition as a func¬ 
tion of temperature for large A has been found. In 
Ref. 58, an unbiased linear response DMFT approach^°^ 
was used to probe stability of the normal phase. Two 
types of instabilities were found, an instability towards 
the solid and an instability towards the excitonic con¬ 
densate. This situation resembles EFKM. Indeed, the 
physics behind formation of the solid and the super¬ 
fluid phases in the strong-coupling limits of EFKM and 
2BHM is similar. As in EFKM, the excitonic instability 
in 2BHM extends to the weak-coupling limit. This is not 
so for the solid phase. In the weak-coupling limit, the 
instability towards solid either does not exist at all or is 
weaker than the antiferromagnetic instability. 

The basic physical properties of the excitonic phase 
were studied in Ref. 102. In Fig. 14, we show a typical 
T-dependence of the magnitude of the order parameter 
\(j)\. The systems selected a polar EC state, (f>iA4>i = 0,^°^ 
which agrees with the antiferromagnetic nn exchange ex¬ 
pected in the the strong-coupling limit (15). The shape 
of |^(T)| is consistent with the mean-field (1 — T/Tc)^^^ 
dependence, expected for DMFT method. As in the 
weak-coupling limit, the excitonic condensation leads to 
opening of a charge gap and appearance of Hebel-Slichter 
peaks in the one-particle spectra, see Fig. 15. This be¬ 
haviour reflects the formal analogy to an s-wave super¬ 
conductor discussed in Sec. II. 

This analogy however, does not extend to the electro¬ 
magnetic properties of the condensate. The neutral exci- 
ton condensate does not contribute to the charge trans¬ 
port, which is facilitated by the quasi-particle excitations. 
In Fig. 15 we show the optical conductivity and dc resis¬ 
tivity at various temperatures. The opening of the charge 
gap leads to an optical gap and exponential increase of 
the resistivity below Tc. 



FIG. 15. Left: The evolution of the one-particle spectral 
functions (diagonal elements) with temperature along the line 
marked in Fig. 14a. The spectra are taken at r=1160, 968, 
921, 892, 829, 725, 580, and 290 K. Right: The corresponding 
optical conductivity. The Inset shows the T-dependence of 
the dc resistivity. Adapted from Ref. 102. Copyrighted by 
the American Physical Society. 

Finally, we discuss the spin susceptibility xV^ i shown 
in Fig. 16. While the normal-phase susceptibility ex¬ 
hibits Curie-Weiss behaviour, the susceptibility of the 
polar EC phase appears T-independent. This observa¬ 
tion holds within the numerical accuracy for the uniform 
susceptibility and approximately also for the local sus¬ 
ceptibility. This behaviour can be understood from a 
single atom picture. In the normal phase, the system is 
locally in a statistical mixture of the LS ground state and 
thermally populated HS multiplet. This leads to Curie- 
Weiss susceptibility and a vanishing of the spin gap. In 
the polar EC phase, the Weiss field mixes the LS and 
HS states and opens a gap of the order U'4>/2 T 
between the local ground state and the excited states. 
The local ground state is a superposition of the form 
a|0) -b /I(| — 1) -b |1)). The spin susceptibility has T- 
independent van Vleck character and a finite spin gap 
appears, Fig. 16. 

Kaneko et used variational cluster approximation 
to study the spin-triplet EC in 2BHM without Hund’s 
coupling (J, J' = 0) and with symmetric bands ta = tb 
for a broad range of interaction parameters. Similar to 
the strong-coupling phase diagram Fig. 10 the authors of 
Ref. 106 found continuous transition between the EC and 
normal state (band insulator), and a first-order transition 
between the EC and AFM phases as shown in Fig. 17. In 
Ref. 55, Kaneko and Ohta extended the study to include 
the Hund’s coupling, which confirmed that J < 0 favours 
the spin-singlet charge-density-wave state, while J > 0 
selects the spin-triplet spin-density-wave state^^’^°^. 

2. Doping 

The weak-coupling theory of doped excitonic insula¬ 
tor was developed by Volkov and collaborators^^^^^ for 
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FIG. 16. Left: The static uniform spin susceptibility Xs 
(circles with error bars) along the line marked in Fig. 14a. 
The stars and dotted line indicate XsiT) for the system con¬ 
strained to the normal state. The red vertical line marks the 
EC transition. Adapted from Ref. 102. Copyrighted by the 
American Physical Society. Right: The imaginary part of the 
local dynamic susceptibility Xioc(c^) at various temperatures 
across the EC transition^®®. 



FIG. 17. The T = 0 phase diagram of 2BHM on a square 
lattice for ta = U = 1 (IF = 4), Vab,Vba, J, J' = 0, A = 6 
obtained with variational cluster approximation. The tran¬ 
sition between the band insulator (N) and the EC phase is 
continuous, the EC/AEM phase transition is of the first or¬ 
der. The thick black line for U' = 0 marks a paramagnetic 
metal. Adapted with permission of authors from Ref. 106. 
Copyrighted by the American Physical Society. 

systems without Hund’s coupling. They showed that 
such systems tend to develop ferromagnetic order due 
to simultaneous appearance of the spin-singlet and spin- 
triplet order. In Sec. HID, we have reviewed the applica¬ 
tion of the weak-coupling approach by Bascones et 
Besides finding Volkov’s ferromagnetic phase (COL) the 
authors of Ref. 85 found also a spin-triplet phase (NC) 
induced by an external magnetic field. The possibility 
of the triplet ferromagnetism was discussed also by Ba- 
lents®® who considered the effect of doping in a strong 
coupling limit on a qualitative level. 

Kunes®® used DMFT to study the effect of doping in 
2BHM with strong Hund’s coupling, which limits the 
possible EC order to spin triplet. The phase diagram, 
shown in Fig. 18, contains the normal phase (open cir¬ 
cles) and three excitonic phases. The polar phase (red), 
discussed in preceding section at half-filling, extends to 




FIG. 18. Phase diagram of 2BHM with hole-doping (counted 
from half-filled case) (left) and chemical potential (right) 
on the horizontal axis. Solid lines mark continuous tran¬ 
sitions, while dashes lines mark the first order transitions 
(right) or boundaries to the phase separation regime (left). 
The coloured dots represent the actual numerical calculations. 
The red colour marks the polar EC phase, the blue and green 
colour mark the FMEC phase. The blue and green phases 
are distinct when the cross-hopping in (7) is strictly zero. 
Adapted from Ref. 65. Copyrighted by the American Physi¬ 
cal Society. 


finite doping levels. It is distinguished from the other 
two phases by absence of ordered spin moment. The re¬ 
lation <f)* A cf> — 0 implies that 4> can be factorised to a 
real vector and a phase factor. The arbitrariness of the 
phase factor is connected to the absence of the cross- 
and pair-hopping in the studied model. The polar phase 
is equivalent to the El phase in Fig. 13. 

At higher doping levels and lower temperatures the sys¬ 
tem enters the FMEC phase with finite magnetisation 
m ~ i{4>* A cf)) 0. The transition between polar and 
ferromagnetic phases proceeds via an intermediate phase 
(blue), which is distinguished from the EMEC only in 
absence of cross- and pair-hopping®®. The EMEC phase 
is equivalent to the COL phase of Ref. 85 (Fig. 13). Sim¬ 
ilar to the T — 0 phase diagram of Ref. 85, first-order 
transitions and phase separation is found at low temper¬ 
atures. It is not clear from the numerical data whether 
a T = 0 polar phase exists at finite doping. At the mo¬ 
ment we can only speculate that the antiferromagnetic 
nn coupling provides a means to stabilise it. A peculiar 
feature of the FMEC is the T-dependence of the magneti¬ 
sation in the vicinity of the continuous transition to the 
normal state, which follows the linear 1 — T/T^ depen¬ 
dence. This behaviour is a consequence of the quadratic 
dependence of magnetisation on the EC order parameter 
m ~ A 4>). 

Finally, we briefly discuss the effect of finite cross- 
and/or pair-hopping, which break the charge conserva¬ 
tion per orbital flavour. As a result the phase factor in 
the polar EC phase is fixed and the system selects either 
the spin-density-wave or spin-charge-density-wave order. 
It is also possible that both types of order are realised 
in parts of the phase diagram. Another consequence of 
finite cross-/pair- hopping will be the absence of a con¬ 
tinuous normal to FMEC transition. The FMEC state is 
selected by terms of the order in the Ginzburg-Landau 
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functional and a continuous transition is only possible if 
the second order terms do not depend on the phase of </>. 
Finite cross-/pair- hopping removes this degeneracy. For 
weak cross hopping one can expect a small wedge of polar 
EC phase separating the normal and FMEC phases. The 
first-order normal/FMEC transitions are still possible. 


IV. MULTI-ORBITAL HUBBARD MODEL 

Little has been done in generalisation of the physics of 
Sec. Ill to systems with more than two orbitals per atom. 
Kunes and Augustinsky^°^ used Hatree-Fock (LDA-I-U) 
approach to study EC in quasi-cubic perovskite with 
nominally 6 electrons in the d shell. Let us demonstrate 
the new features on 5-orbital model describing the d- 
orbital atoms on a cubic lattice. The basic setting of 
such system is similar to half-filled two-orbital model. 
The crystal-field splits the d orbitals into threefold de¬ 
generate t 2 g and twofold degenerate Cg states. With 6 
electrons per atom the lower t 2 g levels are filled and the 
upper Cg levels are empty. An exciton is formed by mov¬ 
ing an electron from a t 2 g orbital to an eg orbital on 
the same atom. Unlike 2BHM with only one possible or¬ 
bital structure of an exciton, in the d-atom there are six 
possible orbital combinations. The cubic symmetry dis¬ 
tinguishes the orbital symmetries of an exciton into two 
3-dimensional irreducible representations Tig and T 2 g. 
Considering the geometry of these two excitons^®^ one 
can infer that Tig excitons are more tightly bound and 
substantially more mobile that the T 2 g ones. Therefore 
only the Tig excitons need to be considered as candidates 
for condensation. The three-fold orbital degeneracy of 
the Tig-excitons adds to the A = I spin degeneracy mak¬ 
ing the order parameter a more complex object. It can be 
arranged to a tensorial form of a 3 x 3 matrix where 
the a indexes the spin and j3 the orbital components. 
thus transforms as a vector under the SO{S) rotations 
in the spin space and as a pseudo vector (Tig representa¬ 
tion) under the discrete symmetry operations of the cubic 
point group. The explicit expression for in terms of 
the on-site occupation matrix can be found in the supple¬ 
mentary material of Ref. 102. The structure of the order 
parameter resembles of superfluid and one can 

use the group theoretical methods applied there to anal¬ 
yse the possible EC phases^°®. More complex structure 
of the order parameter allows for existence of numerous 
distinct phases. 

Another important difference to 2BHM concerns the 
hard-core constraint imposed in the excitons in the strong 
coupling limit. While in the case of 2BHM there cannot 
be more than one exciton on a given atom, in the 5- 
orbital model the hard-core constraint is less restrictive. 
There cannot be more than one exciton of a given or¬ 
bital flavour on atom, but it is possible that two excitons 
with different orbital flavours meet. In the cobaltites 
terminology a single exciton on atom represents the in¬ 
termediate spin S' = 1 state, while the S = 2 high-spin 


state can be viewed as a bi-exciton. The fact that the 
site energy of the high-spin state is lower than of the 
intermediate-spin state^^*^ is then expressed as an attrac¬ 
tion between excitons of different orbital flavour and par¬ 
allel spins. Bosonic models with infinite intra-species, but 
finite inter-species interaction have been studied for two 
species (spinless) bosons®^’®®’"^^’^^ and shown to exhibit 
phases which are not allowed with inter-species hard-core 
constraint. 


V. MATERIALS 
A. Bulk materials 

There have been numerous proposals of materials 
to exhibit excitonic condensation, but very few reali¬ 
sations of EC were actually documented. The early 
candidates for excitonic condensation, which followed 
the weak-coupling picture of proximity to semimetal- 
semiconductor transition, included the group V semimet¬ 
als (Bi, Sb, As) and divalent metals (Ca, Sr, Yb), possibly 
under pressure®®. However, the search for signatures of 
EC in these materials was not successful. 

Wachter and collaborators^^^“^^^ reported observation 
of EC in TmSeo,45Teo.55 under pressure. While they did 
not see a sharp thermodynamic transition, they argued 
that the observed anomalies are consistent with a ’tran¬ 
sition’ in weak (pairing) field. Wakisaka et inter¬ 

preted their photoemission data on Ta2NiSe5 in terms of 

exciton condensation^®®di6.ii7_ Also the lattice distor¬ 
tion and photoemission spectra of a layer compound IT- 
TiSe2 have been interpreted in terms of charge-density- 
wave type exciton condensation and described with the 
weak-coupling theory^^®^^®^. 

While the previous examples involved cases of spin- 
singlet condensation, in the following we will discuss ex¬ 
citonic magnetism. In 1999 Young et reported ob¬ 
servation of 600 K ferromagnetism in Lao.oosCao.ggsBg. 
Ferromagnetism in a slightly doped semiconductor lead 
several groups to generalise the weak-coupling theory of 
excitonic ferromagnetism pioneered by Volkov and col- 
laborators^®^^'* to the case of multiple Fermi surface 
sheets®®’®'*’^®®^^®®. Subsequent investigations showed, 
nevertheless, that the band gap in parent compound 
CaBg is « 1 eV®®® and thus inconsistent with the EC 
scenario. It is now generally accepted that the ferromag¬ 
netism in Laa;Cai_a;Bg arises from defects rather that 
doping. 

Another group of materials where the EC concept 
found its use are iron pnictides®®®"’®®®. The physics of 
these materials is governed by nesting between several 
Fermi surface sheets formed by bands of different or¬ 
bital characters®®®. Several groups studied a simpli¬ 
fied two-orbital model®®’®®®“®®® and its multi-orbital ex¬ 
tensions®®®’®®"® using weak-coupling approaches, and ob¬ 
served a spin-density-wave order with a periodicity given 
by the nesting vector. Arising from nesting between 



19 


bands that mix several orbital characters, the corre¬ 
sponding Weiss field in general couples all possible or¬ 
bital combinations. However, since the nested patches of 
the Femi surface have different dominant orbital charac¬ 
ters a large part of the condensation energy comes from 
orbital off-diagonal pairing, which produces local mag¬ 
netic multipoles but no local moments (polar EC state). 
The generally present, but small, orbital diagonal con¬ 
tributions than give rise to the apparently small ordered 
moments. A first principles calculation of ordered state 
supporting this picture was done by Crincchio et al}^^. 

Kunes and Augustinsky^® proposed that a transition 
observed in some materials of Pra;Cai_2;Co03 (PCCO) 
family^^®^^^® can be understood as an excitonic conden¬ 
sation. Materials from this family exhibit a phase tran¬ 
sition with Tc as high as 130 K which is characterised 
by sharp peak in the specific heat, transition from high- 
T metal to a low-T insulator, disappearance of Co lo¬ 
cal moment response and simultaneous Pr^’*' —)■ Pr'^’*' va¬ 
lence transition. A puzzling feature of the low temper¬ 
ature phase is the breaking of time reversal symmetry 
(in absence of ordered moments) evidenced by Schottky 
anomaly associated with splitting of the Pr^+ Kramer’s 
ground state. The transition to a spin-density-wave EC 
state provides a comprehensive explanation of these ob¬ 
servations and the EC ground state is obtained with 
Hartree-Fock-type LDA-I-U calculations. 

Another class of the materials with potential to ex¬ 
hibit the exciton condensation was proposed by Khali- 
ullin^^°’^^^. He considered a strong-coupling model of a 
d-electron material with cubic crystal field, strong spin- 
orbit coupling and an average d occupancy of four elec¬ 
trons per atom. A scenario possibly realised in materials 
with Re^"*", Ru"'’*", Os"^’*", or Ir®’*' ions. The large crystal 
field restricts the low energy physics to the space spanned 
by t 2 g states. Now the spin-orbit coupling plays the role 
of A in (7) competing with the Hund’s coupling J. Suf¬ 
ficiently strong spin-orbit coupling renders the single-ion 
ground state a singlet S = 0^^^ and the first excited 
state a triplet A = 1. The perturbative treatment of nn 
hopping results in the model similar to (17). There is, 
however, one important difference between (17) and the 
Khaliullin’s model formulated in terms of pseudospin S. 
The real spin S is decoupled from the lattice and thus 
(17) is invariant under the SO{3) spin spin rotations, 
in particular the hopping K± is spin independent. The 
pseudospin S represents a spin-orbital object, which is 
coupled to the lattice and thus the model cannot be in¬ 
variant under continuous pseudospin rotations. This is 
reflected in the hopping amplitudes being S'-dependent, 
i.e., the hopping containing terms known from Kitaev 
modeb^^. 

Finally, the possibility of exciton condensation in lay¬ 
ered cuprates^"^"^ and oxide heterostructures^"*® was dis¬ 
cussed by several authors, but has not been realised so 
far. 


B. Bi-layer structures 

A major direction in the research of exciton condensa¬ 
tion are bi-layer structures'"^®. The basic idea is that a bi¬ 
layer with negligible inter-layer tunnelling provides a sys¬ 
tem where two orbital flavours are to high accuracy inde¬ 
pendently conserved. Exciton condensation takes place 
under suitable conditions including small inter-layer dis¬ 
tance, so that inter-layer electron-electron interaction is 
sufficiently strong, large enough intra-layer electron mo¬ 
bility, so that sufficiently high Tc can be achieved, and 
a matching doping, such that the electron concentration 
in one layer closely matches the hole concentration in 
the other layer. Two types of structures have been stud¬ 
ied. 2D quantum well systems in perpendicular magnetic 
field. The idea here is to use the formation of Landau lev¬ 
els and, in particular, the dependence of the number of 
quantum states per Landau level on the magnetic field as 
a means to achieve the desired electron and holes concen¬ 
trations. The exciton condensation in these systems was 
evidenced by enhanced inter-layer tunnelling^^^ or van¬ 
ishing Hall conductivity^^®. For a review of experimental 
challenges probing exciton condensation in bi-layers we 
refer the reader to specialised literature^"*®. 

Another bi-layer system that has been intensely stud¬ 
ied is the system of two graphene layers separated with a 
dielectric where EC was predicted to take place at room 
temperature*®®. The important and much debated issue 
here is the screening of the inter-layer interaction*®*. The 
system therefore cannot be described with a simple Hub¬ 
bard type model for the strongly correlated electrons. For 
a review of the physics of electron-electron interaction in 
graphene structures we refer the reader to a specialised 
literature*®®. 


VI. OUTLOOK 

Analogy to superconductivity has been the traditional 
driving force of the field of exciton condensation. In par¬ 
ticular, there was a considerable effort in realisation of 
supercurrents driven by the gradient of the condensate 
phase - a hallmark of superfluidity. There are two im¬ 
portant conditions to realise a supercurrent. First, the 
phase invariance (conduction/valence charge conserva¬ 
tion) in the normal phase. Second, the ability to con¬ 
tact separately the electron and the hole parts of the 
exciton. These conditions can be met in bi-layer sys¬ 
tems, where the electrons and hole are spatially sepa¬ 
rated while maintaining sufficiently strong interaction. 
We find it, however, unlikely that similar situation can 
be realised in bulk materials. Already the first condition 
is bound to be violated by various one- and two-particle 
terms in the Hamiltonian as discussed previously. In¬ 
stead of a superfluid with an arbitrary phase the system 
selects either charge(spin)-density-wave or charge(spin)- 
current-density-wave state** with a fixed phase. Omit¬ 
ting possible special points in the phase diagram where 
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the symmetry may be enhanced, e.g., close to a second- 
order boundary between these phases the system will not 
exhibit a superfluid behaviour. 

The main value of the concept of exciton condensa¬ 
tion in bulk materials lies, in our opinion, in providing a 
comprehensive picture of potentially complex phase dia¬ 
grams and understanding long-range orders that are not 
easy to detect because they do not lead to charge(spin)- 
density modulations on inter-atomic scale. Moreover, the 
geometrical form of the order parameter in EC system 
may non be as intuitive as for example magnetisation 
is systems with local moments. For example, exciton 
condensation in s-p systems leads to formation of elec¬ 
tric dipoles® or condensation in a system close to the 
spin-state transition gives rise to a magnetic multipole 
order®®. There are other examples of orbital-off-diagonal 
orders, e.g. orbital currents in cuprates^®®^^®®, nematic 
order in iron pnictides^®®^^®®, hidden order in URu 2 Si 2 ^®® 
where we can only speculate that they may be viewed as 
a special example of the EC physics. It should be clear 
from this discussion that the concept of exciton conden¬ 
sation extended beyond the strict superfluid state is not 
sharply defined and, in particular, the weak-coupling bor¬ 
derline between exciton condensation and ’just a’ Fermi 
surface instability is quite fuzzy. 

We see numerous open questions and possible direc¬ 


tions of further investigation. In the two-band model, 
realisation of spin-current-density-wave phase and in gen¬ 
eral the phase diagram in the presence of cross-hopping 
are to be explored. The Weiss field in the spin-current- 
density-wave state can be viewed as a spontaneous spin- 
orbit coupling - a field that breaks the spin-rotational 
symmetry, but preserves the time reversal symmetry. In¬ 
teraction between the exciton condensate and the lattice 
has attracted attention recently^®®. The multi-band sys¬ 
tems with multiple orbital flavours of excitons are com¬ 
pletely unexplored to our knowledge. A model that can 
shed some new light on the long standing problem of per- 
ovskite cobaltites. 
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